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The differential cross section and analyzing power A y of the pp— >pp7r° reaction have been mea- 
sured at RCNP in coplanar geometry at a beam energy of 390 MeV and the dependence on both 
the pion emission angle and the relative momentum of the final protons have been extracted. The 
angular variation of A y for the large values of the relative momentum studied here shows that this 
is primarily an effect of the interference of pion s- and p- waves and this interference can also explain 
the momentum dependence. Within the framework of a very simple model, these results would 
suggest that the pion-production operator has a significant long-range component. 

PACS numbers: 24.70.+S, 13.75.Cs, 13.60.Lo 



I. INTRODUCTION 

Pion production is the first and probably the simplest 
of the inelastic processes in nucleon-nucleon collisions and 
its understanding may provide us with valuable informa- 
tion about low and medium energy strong interaction 
physics. Motivated by this, over the last decade there 
has been a series of detailed studies of the pp-^ppn re- 
action from near the production threshold up to a proton 
beam energy of 425 MeV. These have been carried out at 
several laboratories, viz. IUCF TSL [3, HQ, and 

COSY-TOF @,0]. 

Initial theoretical calculations 0] for 7r°-production es- 
timated a total cross section that was five times smaller 
than experiment This indicated clearly that some 
essential mechanism was missing from the theory of s- 
wave pion production leading to the S'-state of the fi- 
nal protons. Lee and Riska [3] suggested that this 
might be connected with short-range effects between 
nucleons, and quantitative support for this was found 
in a model with the exchange of a heavy meson cou- 
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pled to the antinucleon-nucleon pair [lOj. On the other 
hand it has been claimed that s-wave pion rescatter- 
ing is not small 0, [l2j and that one can reproduce 
the cross section data without invoking short-range ef- 
fects Studies based on Chiral Perturbation Theory 
(xPTHllBEllIoI show that the contribution of pion 
rescattering is indeed sizable, but that the sign of this 
term is opposite to that of Ref . [TJ [l2j , and this leads 
to an even more severe discrepancy between theory and 
experiment. The convergence of the chiral expansion can 
be seriously questioned for s-wave pion production be- 
cause of the large momentum transfer between nucleons, 
whereas the expansion seems to show convergence in the 
case of p-wave pio n (l7| . An ordering scheme has since 
been discussed [18j . Thus the origin of s-wave pion pro- 
duction is still not clear even as to whether the produc- 
tion mechanism is or is not dominated by short-range 
effects. 

Differential cross sections and polarization observ- 
ables have been measured at several bombarding ener- 
gies [2, |4j, [5j, [Z| and theoretical estimates of these observ- 
able have been made by groups at Jiilich [13] and Os- 
aka [23]. Both approaches include higher partial- waves 
and provide good fits to the total cross section close to 
threshold. However the differential and polarization ob- 
servables are not well reproduced by either group. In 
order to try to identify the origin of the problem more 
clearly it may be useful to attempt a partial-wave de- 
composition. In Ref. , the angular dependences of the 
different polarization observables were developed using a 
general formalism consisting of a complete set of func- 
tions, with coefficients extracted b y fi tting data, and this 
has been extended in later work [2lj . However, due to 
the limited accuracy of polarization data existing at the 
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time, the analysis was done by assuming that only a small 
number of partial- wave contributed and that the momen- 
tum dependence of these amplitude stemmed purely from 
the centrifugal barrier. This latter assumption is very 
doubtful since, for the P-state of final protons, it does 
not lead to a good description of the differential cross 
section Q. More precise data on the angular and mo- 
mentum dependence of the polarization observables, as 
well as of the differential cross section, are highly de- 
sirable. The accumulation of such data may allow one 
to perform a fuller partial- wave analysis of the pp—>ppTr° 
reaction, which would provide greater insight for the the- 
oretical models. 

We report here a measurement of the J)p—>ppir° re- 
action at a beam energy of 390 MeV, where the beam 
polarization is perpendicular to the plane containing the 
detectors. The differential cross section and analyzing 
power, A y , are obtained as functions of the pion emission 
angle, 9 q , in the center-of-mass system and the relative 
momentum, p, of the final protons. The use of a high 
intensity polarized beam in conjunction with a liquid hy- 
drogen target makes it possible to extract the angular 
dependence of A y for different ranges in p. Such data, 
which have not been available previously, might offer se- 
rious constraints on theoretical models. 

If only a few partial- waves are important, the angu- 
lar dependence of the A y x spin-averaged cross section 
is composed of terms proportional to svn9 q and sin20 g . 
The strength of the sin 9 q term is governed by the inter- 
ference of amplitudes corresponding to the production of 
Ps and Pp final states. Here, using the standard nota- 
tion, the final states are labeled by LI, where L and I are 
the angular momenta of the proton pair and the pion, 
respectively. On the other hand, the sin20 g strength is 
determined by S- and higher states. 

The variation with the momentum p, obtained after 
integration over angles, is sensitive to the Ps and Pp pro- 
duction amplitudes. This is studied in a model where it is 
assumed that the momentum dependence of the proton- 
proton wave functions in the final state is the main influ- 
ence for the P-state amplitudes. Furthermore, in a very 
simplistic approach, the wave functions are evaluated at 
a fixed proton-proton separation distance. 

The experimental facility, and in particular the detec- 
tors, polarized beam, and target, are described in Sec.HIl 
Section [HI] is devoted to the steps needed to identify 
and measure the pp-^ppir° reaction. The extraction of 
the observables from data taken under our specific kine- 
matic conditions, where the detection system is essen- 
tially coplanar, is the subject of Sec. IIV1 with the ex- 
perimental results being shown in Sec. |Vj The general 
features of the angular and momentum dependence of the 
cross section and analyzing powers are given in Sec. IVII 
It is shown there that, when comparing our results with 
published data which have large acceptance, it is crucial 
to take account of our particular coplanar geometry. Sec- 
tion IVIII presents the simple phenomenological descrip- 
tion of pion production. Our conclusions are given in 



Sec. lyim 



II. EXPERIMENT 

The experiment was carried out using a 390 MeV po- 
larized proton beam extracted from the cyclotron com- 
plex at the Research Center for Nuclear Physics (RCNP), 
Osaka. The polarized beam was produced by an atomic- 
beam-type polarized ion source with an Electron Cy- 
clotron Resonance (ECR) ionizer [22|, where the polar- 
ization state ('up' or 'down') of the primary beam was 
reversed with a frequency of 1 Hz. The protons were first 
accelerated in the injector AVF cyclotron before being 
further accelerated up to 390 MeV in the main cyclotron 
ring. After extraction, the vertically polarized beam was 
transported to the scattering chamber in the experimen- 
tal hall. Figure Q] shows the top view of the scatter- 
ing chamber as well as of our detection system. The 
scattering chamber containing a liquid hydrogen target 
consisted of an evacuated vertical cylinder. The window 
made of a aramid foil had a horizontal opening angle from 
15° to 110° on either side of the the beam direction. Af- 
ter the target the primary beam was transported to the 
beam dump where a Faraday cup monitored the beam 
intensity for the two polarization states. The beam cur- 
rent was limited to around 1 nA in order to minimize the 
dead-time of the data acquisition system. Under these 
conditions the system, the details of which are described 
in Ref. [23j, had an efficiency of between 75% and 87%. 

A. Detector 

The detector system was symmetric with respect to 
the plane containing the beam axis and the direction of 
the beam polarization. The two outgoing protons from 
the pp^ppir reaction were detected simultaneously in 
an array of plastic scintillators. The measurement of two 
pairs of polar and azimuthal angles and two kinetic en- 
ergies is sufficient to identify the pp—>ppir° reaction and 
determine the five independent kinematic variables. The 
system covered laboratory polar angles 15° < 9 < 35°. 
The maximum angle of protons from the pp^>ppir° re- 
action at 390 MeV (32.1°) is well inside the angular ac- 
ceptance. The minimum angle was limited for both pro- 
tons such that only proton pairs with relative momenta 
from 150MeV/c up to the kinematic limit of 220MeV/c 
could be registered with this system. The detector was 
therefore well suited for the investigation of the high rel- 
ative momentum region, where the pion-production am- 
plitudes leading to P-state proton pairs should have their 
maximum strength. However, it is important to note that 
the cm. polar angle of the relative momentum vector for 
final protons was confined to 90° ± 30° with respect to 
the beam axis. 

The energy of a scattered proton was deduced from the 
amplitude signal in one of the P-counters. Such a counter 
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FIG. 1: Top view of the layout of the experiment. The shaded region shows the walls of the target chamber and the dashed 
lines its windows. In the right-hand panel a three-dimensional sketch shows the positions of the hodoscope and _E-counters. 



consists of a set of five plastic scintillators in which each 
element has a trapezoidal shape, with front and back 
faces of area 35 x 35 mm 2 and 60 x 60 mm 2 respectively. 
The 350 mm length is sufficient to stop protons with ener- 
gies up to 250 MeV so that all protons from the pp^ppir 
reaction at 390 MeV were stopped. The front face of 
each counter was positioned at 500 mm from the target. 
The scintillator hodoscopes in front of the counters were 
used to determine the directions of outgoing particles; a 
single 3 mm thick element of area 17.5 x 30 cm 2 covered 
±17mrad horizontally and ±30mrad vertically Two el- 
ements were placed in front of each of four -©-counters 
with a 1 mm overlap, though only one element was used 
for the counter at the largest angle. There was thus a to- 
tal of nine elements on each side of the beam. The angle 
of the detected particle was denned by the center of each 
element. 

The anticoincidence counter ( marked 'Anticounter' in 
Fig. [1} was used to eliminate accidental events coming 
from elastic proton-proton scattering. This covered the 
angular range 45° — 72°, where recoil protons from the 
elastic scattering hit the counter in combination with 
the scintillator hodoscopes. The information from the 
counter was used in the offline analysis and about 90% 
of the accidentals could be suppressed in this way. 

One hodoscope element and a scintillation counter 
placed at 60°±1° (denoted as the 'Polarimeter' in Fig. |lj) 
were used to monitor the beam polarization. The fast 
and recoil protons from elastic proton-proton scattering 
were detected in coincidence. According to the SAID 
database (zjj . the pp^pp analyzing power at this en- 
ergy and angle is A y — —0.364 ± 0.007, where the er- 
ror bar has been obtained by looking at typical data in 
this region. The average values of the beam polariza- 
tions deduced on this basis were P-\ — 0.66 ± 0.01 and 
Pj = 0.70 ± 0.01, where the arrows indicate the spin 
states of proton beam and the errors are statistical. The 
polarizations were quite stable during the experiment and 
all values lay within about ±0.03 of the averages. 



B. Target 

The liquid hydrogen (LH2) target was 8.5 mm thick, 
with windows made of 25 /jm thick aramid foil [25| . Its 
temperature was controlled in the range of 14 K to 20 K 
and kept stable at 17.6 ± 0.4 K during experimental runs. 
The empty target runs, which were carried out with hy- 
drogen gas at temperatures between 22 K and 27 K, were 
used to estimate the background from other construction 
materials as well as from the residual gas frozen on the 
windows. 



C. Data taking 

A signal was registered when a charged particle hit 
both of the scintillator hodoscopes and the correspond- 
ing ©-counter. For the measurement of the pp—>ppir° 
reaction, the trigger conditions were set so that both 
right- and left-side counters were required to give sig- 
nals in coincidence. In addition, measurements without 
the coincidence requirement were also performed in or- 
der to detect single protons from elastic proton-proton 
scattering. Such measurements were also performed us- 
ing an unpolarized proton beam in order to check the 
determination of the luminosity as well as of the beam 
polarization. These are the subject of the following sec- 
tion. 



III. DATA ANALYSIS 

A. Particle identification 

Particle identification was achieved by the AE-E 
method. The amplitude signal, related to the energy loss 
in the scintillator hodoscope, is plotted in Fig. [5] against 
the amplitude signal from the ©-counter. The heavily 
populated band arises from protons, associated mostly 
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with 7T° production, that stopped inside the counter. The 
pp^ppir events of interest are well separated from the 
two lower islands. These are produced by elastically scat- 
tered protons, which have energies above 265 MeV and 
thus do not stop inside the counter, and positive pions 
generated through the pp^pnir^ reaction. 
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FIG. 2: Identification of protons from the pp— >ppn° reaction. 
The signal from the hodoscope is plotted versus that from 
the _E-counter. The intense band is due to protons that have 
stopped in the scintillator. Groups of events corresponding 
to elastically scattered protons and pions from the pp—>pnn + 
reaction are indicated and separated by the cut shown as the 
solid line. 
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FIG. 3: Two-dimensional plot of the proton energy T2 de- 
posited in the right-hand .E-counter versus the same quantity 
Ti for the left-hand counter. The hodoscope elements, lo- 
cated closest to the beam line on both sides, were selected. 
Panels (a) and (b) show the simulated and experimental data, 
respectively. 



B. Energy calibration of the i5-counter 

At a beam energy of 390 MeV, the energies of the de- 
tected protons from the pp—>ppir° reaction varied be- 
tween 40 MeV and 226 MeV. The energy dependence of 
the amplitude signal from the ^-counter could be cali- 
brated above 100 MeV using elastic proton-proton scat- 
tering events. The hodoscope and i?-counter were set 
to cover the angular range 40° — 60° and, for an angle 
fixed by the hodoscope, the monoenergetic recoil proton 
was measured by an E-counter element. The energy res- 
olution was found to be better than 2% (FWHM) for a 
proton energy of 200 MeV. 

Below 100 MeV, information from the pp—>ppir° reac- 
tion itself was used. Figures^ and[3}} show, respectively, 
the simulated and measured energy correlation between 
protons that hit a pair of hodoscope elements to the right 
and left of the beam. The locus of ir° event can be clearly 
seen and this was used to extract a data sample where 
the energy in the left (right) counter was higher than 
100 MeV, while that in the right (left) was lower, as in- 
dicated by Case 1(2) in Fig. (3^- In addition a pedestal 
value from the ADC module was also used at an energy 



of 17 MeV, corresponding to the minimum energy of pro- 
tons detected in the -E-counter. In this way, the am- 
plitude signals for all combinations of right and left-side 
counters were calibrated using third-order polynomials 
with an accuracy of better than 4%, as judged from the 
position and width of the resulting ir° missing-mass peak. 



C. Identification of the pp^pp7r° reaction 

Having measured the energies and angles of the two 
protons, good events were selected on the basis of the 
missing-mass spectrum. A clear ir° peak, with a width 
of 7.4MeV/c 2 (FWHM), is seen in Fig. H which shows 
the totality of events obtained with the LH2 target. The 
background from the target foil was subtracted by utiliz- 
ing the empty target runs normalized to the integrated 
beam intensity. The contribution from the random coin- 
cidences was estimated by considering outgoing protons 
from different beam bunches. The correspondence of the 
peak position with the mass of the ir° is consistent with 
the energy calibration of the counters. 
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FIG. 4: The pp—>ppX° missing-mass spectrum obtained for 
all geometric configurations using the LH2 (solid histogram) 
and the gaseous hydrogen target of the 'empty run' (dotted 
line). The estimate of random coincidence events is shown by 
the dashed line. The peak corresponding to 7r° production has 
a mass of 135.0 MeV/c 2 and a width of 7.4 MeV/c 2 (FWHM). 



For all pp^ppn candidates, the pion emission angle 
9 q and the relative momentum p and its angle 6 P in the 
overall center-of-mass system, were reconstructed on an 
event-by-event basis. The data were then grouped in nine 
intervals in 9 q and three in p, with a finer binning in p 
being used to study the momentum dependence of the 
analyzing power. 

Missing-mass spectra were constructed for each com- 
bination of bins and spin states and used to extract the 
yield of tt° events. The backgrounds from random coin- 
cidences and target foils were subtracted before showing 
the data in Fig. [5] for typical conditions. There is some 
residual background for 6 q < 20°, whose shape varies 
with the kinematic conditions. This arises from 'non-full' 
events where, in contrast to 'full' events, recoil protons 
from 7T° production do not deposit their full energy in 
the -Encounter due to a nuclear reaction in the scintilla- 
tor material. 

The shapes of the distributions for non-full and full 
events were simulated in a Monte Carlo program based 
on GEANT3 [2(1, which took into account the energy 
resolution of the scintillator material, as well as hadronic 
and electromagnetic interactions. The simulated events 
were passed through the same analysis chain as the mea- 
sured data and fitted to the results shown in Fig. 03 The 
normalization of the non-full events was determined from 
the tails of the distributions. In order to check the sys- 
tematic uncertainty in the treatment of the non-full back- 
ground, the 7r°-yield was also determined by selecting all 
events around the pion peak without any subtraction and 
the corresponding uncertainty is discussed in Sec. [VJ 
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FIG. 5: The pp^ppX missing-mass spectra for different se- 
lected kinematic regions, with the relative momentum p being 
indicated at the top of the panels and the angle of the pion 
8 q being shown within each. The experimental data are pre- 
sented after subtracting the backgrounds from the random 
coincidences and target window. The hatched regions show 
normalized distributions of 'non-full' events, where not all of 
the energy is deposited in the scintillator. The solid lines show 
the sums of these plus the good events. The shapes of both 
distributions have been determined through simulations with 
GEANT3 [H 



D. Detection efficiency of i5-counters for 'full' 
events 



The detection efficiency for full events, where protons 
deposit their full energy inside an Encounter, was esti- 
mated using the data obtained for the energy calibration 
of the ^-counters discussed in Sec. IIIIB1 The efficiency 
varies between 0.85 and 0.65 over the measured energy 
range of 100 MeV to 200 MeV. The efficiency obtained by 
the Monte-Carlo simulation was checked by comparing it 
with that obtained from the measurement. Both results 
are in agreement to within 15% and this is included in 
the overall uncertainty quoted for the cross section. 



E. Luminosity 

The absolute value of the luminosity was determined 
from measurements of the beam intensity and the density 
of the LH2 target and this gave an integrated luminosity 
of 16.1 pb _1 . The value was verified by measuring elas- 
tic proton-proton scattering by identifying single protons 
through the AE-E method. The unpolarized pp^pp dif- 
ferential cross sections obtained in this way are shown 
separately in Fig. [5] for the same angular interval to the 
left and right of the beam direction. The results agree on 
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FIG. 6: The laboratory differential cross sections measured 
for elastic proton-proton scattering using hodoscopes placed 
to the left and right of the beam. The lines are predictions 
taken from the SAID database 12411. 
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FIG. 7: The analyzing power for elastic proton-proton scat- 
tering measured by detecting single protons with hodoscopes 
placed to the left and right of the beam, using beam polar- 
izations determined by the coincidence method. The lines are 
estimates taken from the SAID database 12411 . 



average to within ±2% with the predictions of the SAID 
program [24| . However, there are fluctuations of up to 
±5% around the predicted curves and these deviations 
were included in the determination of the acceptance as 
corrections associated with the geometrical uncertainties 
of individual hodoscope elements. 



F. Beam polarization 

The beam polarizations were also checked by measur- 
ing elastic proton-proton scattering where the selection 
criteria for a single proton is the same as that discussed 
previously. The analyzing powers obtained in the offline 
analysis on the basis of the already determined beam 
polarizations P^ are shown in Fig. [7] for the same lab- 
oratory angular range to the left and right of the beam 
direction. The results are consistent with those given by 
the SAID program [24[ and the systematic uncertainty 
in the determination of the beam polarizations, includ- 
ing that coming from the SAID database, was found to 
be below 6%. 



IV. OBSERVABLES 



Definition 



Five independent variables £ are required to describe 
the three-body final state and these we take from the rela- 
tive momentum p of the two final protons and that of the 
pion q in the overall cm. system. Since the magnitudes 
of these momenta are linked by energy conservation, the 
resulting set of variables consists of the magnitude of p 
plus two polar angles 9 with respect to the beam axis and 
two azimuthal angles 0, i.e. £ ={8 q ,(p q ,6 p ,tp p ,p}. The 
corresponding differential cross section will be denoted 



da 



dflqdflpdp 



(1) 



In terms of the Cartesian observables, the dependence 
of the cross section on the vertical polarization P y is given 
by 



a(Z,P y ) = a (0 [l + P y A y (0], 



(2) 



where CTo is the spin-averaged cross section. 

Our goal is to obtain the spin-averaged cross section 
and analyzing power at the polar angle of 8 p — 90° in 
coplanar geometry, i. e. ip p — ip g = 0° . These observables 
are obtained from the experimental data through 
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(P^JVl + P^N 



{LjNi — L^N])/ (costpg) 
(PiL^ + P T i T 7V|) 



(3) 
(4) 



where N^i is the spin-dependent yield, Pf^ the beam 
polarization, and [ the luminosity. The detection effi- 
ciency e includes that of the data acquisition system and 
the acceptance of the detector system, as determined by 
the Monte-Carlo simulation using the phase-space model. 
The average value of the cosine of the pion angle, (cos <p q ) , 
is also determined through the Monte-Carlo simulation, 
as discussed in Sec. IVl 



B. Angular dependence 

For the later discussion it is convenient here to de- 
scribe the angular dependence expected for the unpolar- 
ized cross section and analyzing power of the pp^ppir 
reaction. These, as well as other observables, have been 
discussed in terms of partial- wave amplitudes in Ref. 0] , 
where total orbital angular momentum excitations up to 
L pp + = 2 were considered. Taking (p p = 0° and 
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O p — 90°, the relevant formulae reduce to 

a (9 q ,ip q ,p) = E + F 1 +H°°-(H™ + F 2 +K) + 
(I + - H$°) (3 cos 2 9 q - 1) + H%° sin 2 9 q cos 2ip q , 

(5) 

* A v (9 q ,<p q ,p) = {(Gf -Gf + Gf) sin 9 q 

+ (P° + Hf - Hf + Hf) sin 26 q } cos Vq , 

(6) 

where each of the coefficients is a function of the relative 
momentum p. The coefficients E, F, and H represent the 
absolute-squares of amplitudes leading to Ss, Ps, and 
Pp states, respectively, whereas those with G, I, and K 
reflect the interferences of type PsPp, SsSd, and SsDs 
respectively [2(. 

V. RESULTS 

The measured values of the spin-averaged cross section 
and the analyzing power are shown in Table U as func- 
tions of the pion polar angle q for three ranges in the 
relative momentum p. The mean values and uncertain- 
ties in the determination of 6 q have been estimated from 
the Monte Carlo simulations, where possible fluctuations 
in the beam energy and geometrical uncertainties in the 
detector system were taken into account. 

The combinations of up-down beam polarizations and 
left-right detectors provide two measurements of A y as 
well as of the spin-averaged cross section and these should 
be consistent. However, due to the uncertainty of the en- 
ergy calibration of the -E-counters on the two sides, the 
data show systematic differences between these combi- 
nations. The deviations from the mean value were in- 
cluded in the systematic error for an individual angular 
bin. Compared with this, the uncertainty due to the 
treatment of the non-full events as discussed in Sec. If 11 Cl 
is negligible for A y , though it is the dominant systematic 
error for the spin-averaged cross section. 

The size of the 9 p angular bin is quite large (90° ± 30°) 
and so the analysis was repeated using the smaller angu- 
lar interval (90° ±10°) in order to estimate the resulting 
uncertainty. The results for the spin-averaged cross sec- 
tion agree to within 5%, whereas the analyzing power 
results are in agreement to within the statistical errors. 
The overall systematic uncertainty for the spin-averaged 
cross section was estimated to be less than 16% when 
the uncertainties from the efficiency of full events and 
the angular bin were taken into account. As mentioned 
in Sec. IIII F[ the overall systematic uncertainty in A y 
coming from the determination of the beam polarization 
is believed to be less than 6%. 

It is important to note that, due to the finite size of our 
counters, the data were not taken strictly at azimuthal 
angles (p q — 0°(180°) and the effects of this angular 
spread increase when 9 q approaches 0°(180°) and p gets 



close to its maximum allowed value. In Table Q] we show 
estimates of the average values of (cos<p q ) and (cos2</? 9 ) 
evaluated using our simulation. It is seen from Eq. ([5]) 
that the product o§A y is proportional to cos ip q but is 
independent of cos2(/? 9 and so the estimates of (cos<p q ) 
have been used in Eq. (J3]) to deduce A y . There is an 
explicit cos2tp g in Eq. |5]). which means that the value 
of do that we have measured is different from that at 
ip q — 0°. As a consequence, when dividing aoA y by <tq 
to derive the analyzing power A y , some ip q dependence 
still remains. 

To quantify the changes caused by the tp q variation, we 
define a correction factor dt using Eq. ([5]), 

d t = cr ({cos2(p q }) /a (cos2(p q = 1) . (7) 

The d t given in our Table [J were obtained for our kine- 
matic conditions by interpolating the coefficients quoted 
in Table IV of the IUCF work 0. It must be noted that 
the value of H®° appearing in Eq. could not be de- 
termined in the IUCF work and the value is assumed to 
be zero in the analysis. The calculations show that the 
cross section increases as <p q moves away from 0° and 
hence that the analyzing power decreases. 

The spin-averaged cross sections and analyzing powers 
found by integrating over the polar angles 6 q are shown 
as functions of the relative momentum p in Table HT1 The 
mean values and uncertainties in the determination of p 
have been estimated in the same way as for the pion polar 
angle 9 q in Table HI The spin-averaged cross sections and 
analyzing powers have been obtained at average values 
of cos 2(p q which depend upon cos 6 q and they must be 
subjected to a correction factor d t , similar to that given in 
Tablefl] in order to extrapolate the results to cos 2(p q = 1. 

VI. DISCUSSION 

A. Angular dependence 

The variation of the spin-averaged cross section with 
pion angle is shown in Fig. [5] for the three momentum 
ranges. Since, for this purpose, the two initial protons are 
identical, the data are presented as functions of cos 2 6* g . 
The variation seems to be linear, as expected on the basis 
of Eq. {5J and the data were therefore fit with the form 

cr = a (1 + ccos 2 6 q ) . (8) 

The resulting parameters with and without the dt factor 
are given in Table UTTl In all cases the % 2 /d.o.f. were close 
to unity. The results, with and without the dt modifica- 
tion, are both consistent with c being constant. 

The angular distribution of the pion has been inves- 
tigated by several experimental groups in the 400 MeV 
region, though with somewhat conflicting conclusions re- 
garding the slope parameter c. Thus c = 0.28 ± 0.20 
was found at IUCF @, 0.19 ± 0.01 at TSL (PROMICE- 
WASA) [|, -0.35 ± 0.03 at TSL (WASA) @, and 
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TABLE I: The spin-averaged pp^ppir cross section ao and analyzing power A y measured in different ranges of relative 
momentum p at the polar angle of 9 P = 90° with its angular bin of ±30° . The first and second error bars in ao and A y show, 
respectively, the statistical and systematic uncertainties, whereas that for the pion emission angle 9 q is the standard deviation 
that includes the effects of the binning. The overall systematic uncertainty for A y is estimated to be less than 6% and that 
for (To to be less than 16% and these are not included in the systematic uncertainties for an individual angular bin. The 
observables have been obtained at the average values of (cos 2ip q ), whereas the value of (cos</9 9 ) is involved in the evaluation 
of A y . In order to get estimates of the values at cos 2^ = 1, the cross section should be divided by the correction factor dt 
and the analyzing power multiplied by the same factor, as discussed in the text. 



e q 


P 


ao 


Ay 


(COS <fiq) 


(cos 2tp q ) 


d t 


(deg) 


(MeV/c) 


(nbsr~ 2 (MeV/cT 1 ) 










14 ±7 




7.64 ±0.03 ±0.16 


-0.085 ± 0.008 ± 0.005 


0.653 


0.206 


1.01 


31 ±7 




7.26 ±0.03 ±0.19 


-0.191 ±0.006 ±0.018 


0.950 


0.810 


1.01 


50 ±7 




5.84 ±0.02 ±0.23 


-0.259 ± 0.007 ± 0.005 


0.977 


0.910 


1.02 


70 ±7 




4.89 ± 0.03 ± 0.38 


-0.345 ± 0.009 ± 0.047 


0.984 


0.938 


1.02 


90 ±7 


160 - 180 


4.92 ± 0.03 ± 0.23 


-0.362 ±0.018 ±0.066 


0.985 


0.943 


1.02 


110 ±7 




5.24 ±0.03 ±0.13 


-0.307 ±0.011 ±0.018 


0.982 


0.933 


1.02 


130 ±7 




5.95 ±0.03 ±0.09 


-0.270 ± 0.009 ± 0.022 


0.972 


0.896 


1.02 


148 ±7 




7.32 ±0.03 ±0.31 


-0.208 ±0.007 ±0.009 


0.941 


0.779 


1.02 


166 ± 7 




7.57 ±0.03 ±0.10 


-0.084 ±0.009 ±0.005 


0.590 


0.101 


1.01 


16 ±8 




5.56 ±0.02 ±0.18 


-0.076 ±0.011 ±0.024 


0.536 


0.072 


1.01 


32 ±8 




5.26 ± 0.02 ± 0.18 


-0.167 ±0.006 ±0.004 


0.905 


0.679 


1.02 


51 ±8 




4.49 ± 0.02 ± 0.03 


-0.259 ± 0.007 ± 0.005 


0.959 


0.843 


1.02 


71 ±8 




3.74 ± 0.02 ± 0.16 


-0.279 ± 0.010 ± 0.023 


0.972 


0.890 


1.03 


90 ±8 


180 - 200 


3.38 ± 0.02 ± 0.14 


-0.339 ± 0.008 ± 0.024 


0.974 


0.902 


1.03 


109 ±8 




3.69 ± 0.02 ± 0.12 


-0.310 ±0.008 ±0.029 


0.971 


0.886 


1.03 


129 ±8 




4.08 ± 0.02 ± 0.19 


-0.215 ±0.007 ±0.006 


0.956 


0.833 


1.03 


148 ±9 




5.09 ± 0.02 ± 0.12 


-0.150 ±0.006 ±0.007 


0.904 


0.661 


1.02 


163 ±8 




5.76 ±0.02 ±0.02 


-0.096 ± 0.009 ± 0.006 


0.561 


0.043 


1.01 


22 ± 14 




2.34 ±0.02 ±0.06 


-0.114 ±0.031 ±0.042 


0.366 


-0.129 


1.02 


37 ± 12 




2.04 ±0.02 ±0.05 


-0.143 ±0.017 ±0.010 


0.725 


0.303 


1.03 


54 ± 12 




1.77 ±0.02 ±0.02 


-0.194 ±0.015 ±0.011 


0.833 


0.533 


1.04 


72 ± 11 




1.57 ±0.01 ±0.09 


-0.253 ±0.016 ±0.025 


0.881 


0.646 


1.05 


91 ± 11 


200 - 220 


1.48 ± 0.01 ± 0.08 


-0.255 ±0.016 ±0.038 


0.896 


0.677 


1.05 


109 ±11 




1.54 ±0.01 ±0.12 


-0.233 ±0.015 ±0.028 


0.894 


0.664 


1.05 


127 ± 13 




1.67 ±0.01 ±0.13 


-0.168 ±0.015 ±0.006 


0.848 


0.546 


1.04 


143 ± 13 




1.89 ± 0.02 ± 0.10 


-0.164 ±0.016 ±0.011 


0.730 


0.305 


1.03 


157 ± 15 




2.04 ± 0.02 ± 0.22 


-0.120 ±0.032 ±0.066 


0.349 


-0.130 


1.02 



-0.19±0.02 at COSY-TOF |7|. In the last of these exper- 
iments, a dependence of the parameter upon the relative 
momentum could be established, with c changing from 
negative to positive as p increases. However, even the 
TOF result of 0.09 ± 0.01 obtained for p > 160 MeV/c is 
significantly smaller than ours for a similar momentum 
range. This difference can be explained in terms of the 
different geometries of the two experiments. 

Whereas at COSY-TOF the full phase space was ex- 
plored, our data were taken in essentially coplanar condi- 
tions, i.e. ip p — ip q — 0°, and with the proton polar angle 
of 9 p = 90°. Using Eq. with the coefficients taken 
from Table IV of the IUCF work 0] , one can derive a re- 
lation between the slope parameter (c) for full acceptance 
and the value (c*) to be found with our limited coverage. 
Since H®° is negative Q, c* is always larger than c and 
for p > 160 MeV/c one expects that c* « 0.3 ± 3.0c. Us- 
ing this relation with the COSY-TOF value of c, a value 
of c* — 0.57 ±0.03 is predicted and this is consistent with 
our results. 



The sign of the slope parameter c is positive for high 
p, reflecting the importance there of the Pp contribution. 
On the other hand, both the TSL [| and COSY-TOF Q 
data show that c is negative for p < 53MeV/c, proba- 
bly due to a SsSd interference that falls very fast with 
increasing p. 

Figure[9]presents the dependence of A y on the pion po- 
lar angle for the three relative-momentum regions. The 
solid line shows the fits resulting from using the general 
angular dependence of Eq. © , viz. 

gqA v — a sin 9 q + b sin 29 q , (9) 

where the values of the free parameters a and b are listed 
in Table ITlTI In the fitting, the uncertainty in the angular 
determination has been included in the \ 2 minimization. 

It is seen from Table IIIII that b is consistent with zero 
for all three momentum bins. Since the parameters a 
and b correspond, respectively, to the PsPp and (SsSd, 
PpPp) contributions in Eq. ([5]), this shows clearly that, 
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TABLE II: The pp—>ppir spin-averaged cross section, 
da/difiqdripdp, after integration over 8 q , presented in differ- 
ent bins of the relative momentum p. The corresponding 
integrated values of the analyzing power A y are also given. 
The first and second error bars denote the statistical and 
systematic uncertainties, respectively, whereas that on the 
relative momentum p is the standard deviation that includes 
the effects of the binning. The observables were measured 
over a range of ip q and, in order to get estimates of the values 
at cos2<^) 9 = 1, the cross section should be divided by the 
correction factor dt and the analyzing power multiplied by 
the same factor, as discussed for Table [T] 



p da/difiqdttpdp 


Ay 


d t 


(MeV/c) (nb(radsrMeV/c)~ 1 ) 







155 ±4 12.33 ±0.04 ±1.37 

165 ±4 12.39 ±0.03 ±0.28 

175 ±4 10.88 ±0.03 ±0.24 

185 ±4 9.52 ±0.02 ±0.13 

195 ±4 7.42 ±0.02 ±0.19 

205 ±4 4.76 ± 0.02 ±0.15 

214 ±4 2.12 ± 0.01 ±0.13 



-0.278 ± 0.005 ±0.014 1.02 

-0.283 ±0.004 ±0.018 1.02 

-0.266 ±0.004 ±0.011 1.02 

-0.251 ± 0.004 ± 0.005 1.02 

-0.223 ± 0.004 ± 0.010 1.03 

-0.203 ±0.006 ±0.011 1.04 

-0.136 ±0.010 ±0.040 1.03 
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FIG. 9: The dependence of the ~pp^>ppir° analyzing power 
on the pion polar production angle 8 q for the three regions 
of relative momenta of the final protons. The error bars in- 
clude the statistical and systematic uncertainties. The solid 
curves are fits of Eqs. ((8} and ((9} to these and the data of 
Fig. [8] with the resulting parameters being given in Table [TTTl 
The upper and lower dotted curves in the lowest momentum 
region have been calculated from Eqs. (0) and (|6]) using the 
coefficients given in Ref. |2| at 375 and 400 MeV, respectively. 
The curves bordered by shaded areas are the theoretical re- 
sults from the Osaka group, with the shading indicating the 
expected uncertainties in the calculation [2pj ] . 



FIG. 8: The angular dependence of the spin-averaged 
pp—>ppiv cross section for the three regions of relative mo- 
mentum indicated. The error bars include the statistical and 
systematic uncertainties. The solid lines show the results of 
fitting the data with the linear form of Eq. (0, with the re- 
sulting parameters being given in Table [TTTl 

for the large values of p investigated in our experiment, 
the PsPp interference terms dominate. 

If the values of all the coefficients given in Table IV of 
Ref. [|J are inserted into Eqs. (O and ([6]), these predict 
the angular dependence shown in Fig.[9K. The upper and 
lower dotted lines are obtained by using the numbers cor- 
responding to 375 and 400 MeV, respectively. It must be 



noted that the IUCF parameters were obtained from av- 
erages of A y over the full range of allowed proton-proton 
relative momenta Q. The close agreement found in the 
160-180 MeV/c bin when using the 400 MeV coefficients 
is due, in part, to the cross section being maximal in this 
region. 

In contrast, the model calculations of the Osaka 
group [201 ] . shown in Fig. [51 demonstrate a much more 
asymmetrical behavior than the experimental data and 
greatly underestimate the magnitude of the analyzing 
power. The predictions for the Pp or Ps transitions are 
therefore much smaller than those found experimentally. 
More theoretical studies are needed in order to elucidate 
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TABLE III: Parameters extracted by fitting the angular dependence of (Jo(O q ) and ooA v (Q q ) with Eqs. {(SJ) and ©, respectively. 
The values in brackets show the results obtained when the correction factor dt is included in the fitting process. 



V 


a 


c 


a 


b/a 


(MeV/c) 


(nb/(sr 2 MeV/c)) 




(nb/(sr 2 MeV/c)) 




160 - 180 


4.82 ±0.10 (4.70 ±0.10) 


0.62 ±0.04 (0.64 ±0.05) 


-2.00 ±0.12 


0.05 ±0.06 


180 - 200 


3.47 ± 0.05 (3.35 ± 0.06) 


0.72 ±0.03 (0.76 ±0.04) 


-1.23 ±0.06 


0.04 ± 0.07 


200 - 220 


1.44 ± 0.04 (1.37 ±0.04) 


0.64 ±0.08 (0.71 ±0.09) 


-0.40 ±0.03 


0.06 ±0.11 



20 



(a) 



previously remarked, the cross section seems to be max- 
imal at our lowest values of p, as does the magnitude of 
the analyzing power. 
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VII. PHENOMENOLOGICAL ANALYSIS 
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It was already stressed that the combination cr§A v is 
not sensitive to the cos2c/? g term in Eq. §5§ and so does 
not suffer from the resultant d t ambiguity. Since the 
integral of the sin 29 q term in Eq. © over the pion angle 
vanishes, the dependence of the integrated gqA v {p) can 
be expressed in terms of the interference between the two 
types of P-state amplitudes multiplied by the three-body 
phase-space factor p(jp)\ 



a A y (p) = PsPpp{p). 



(10) 



-0.3 



The PsPp corresponds to a linear combination of the 
coefficients Gf° 2 4 in Eq. ©. According to Ref. Q, G\° , 
which does not contain any 8 p dependence, contributes 
about 95% of the total magnitude of A y so that to a good 
approximation the G\^ A may be neglected. In this case, 
PsPp results from s-p interferences with the same P- 
states of the final protons. It can therefore be expressed 
as 



160 180 200 220 

p [MeV/c] 

FIG. 10: (a) Spin-averaged cross section as a function of the 
relative momentum of the final protons, (b) The relative mo- 
mentum dependence of the integrated analyzing power. The 
error bars include the statistical and systematic uncertainties. 

the origin of this disagreement. 

B. Momentum dependence 

Figures [TUb and [TUb show, respectively, the spin- 
averaged cross section and analyzing power integrated 
over 6 q as functions of the relative momentum of the fi- 
nal protons. The error bars include both systematic and 
statistical uncertainties. On kinematic grounds, A y must 
tend to zero as p approaches the highest allowed values 
because the pion momentum vanishes in this limit. As 



PsPp = C /3 PqS f 3 p op + C 2 f3p. 2S /3p 2 p, (11) 

where Co, 2 are parameters and f3 P<J 2S and /3 Po 2P repre- 
sent the pion s- and p-wave amplitudes, respectively. 

A partial-wave amplitude for pion production may be 
approximated in terms of the overlap integral involving 
the pion-production operators 11(7") and the radial wave 
function of initial Ui(r), and final u^(r) protons and the 
pion plane wave ji(^qr), 

ffi= Ju^j^q^fl^u^^dr, (12) 

where r is the proton-proton separation distance, and 
I and q is the pion angular momentum and momen- 
tum, respectively. Analytic forms for the Ss-wave pion- 
production amplitude may be found in Ref. [To| . 

If the pion-production operator is taken to be 
momentum-independent, any strong momentum depen- 
dence of the partial-wave amplitude must be ascribed to 
the variation of the radial wave function of the final pro- 
tons and pion with p and q. To investigate this further, 
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FIG. 11: The relative momentum dependence of the angle- 
integrated A y ao(p). The error bars include the statistical 
and systematic uncertainties. The lines show the variation 
predicted for the 3 p2 state when using Eq. (|13[1 with different 
values of the interaction range. The solid line shows the calcu- 
lation with R — 1 fm, whereas the dashed line and dot-dashed 
line show the calculations with R = 2fm and R = 3fm, re- 
spectively. 

we approximate Eq. (fT2|) by evaluating the product of 
the wave functions at a fixed distance r = R, in which 
case 

f n <xvF f {R)ji{\qR). (13) 

The 3 Po,2 radial wave functions have been derived from 
the Paris potential with the Coulomb interaction [27j |. 
Since the experimental data cannot distinguish between 
the contributions from the 3 Po and 3 p2 states, the calcu- 
lations have been carried out taking them into account 
one at a time. However, there is very little difference be- 
tween the two sets of results and only those for the 3 P2 
states are shown in Fig. 111! The overall normalization 
factor Co, 2 has been fixed on the basis of Eq. ([6]) to sat- 
isfy J a A y (p) dp = crtot {Gf -Gf +Gf)/Wn, where the 
value crtot = 65 /ib has been obtained by interpolation of 
the IUCF data 0. 

The <Jf)Ay data are only sensitive to a value of the dis- 
tance R = Rp s = Rp p . As demonstrated in Fig. [TTJ the 
curves obtained from Eq. ([13")) using radii R = 1, R = 2, 
and R = 3 fm show significantly different dependences on 
the relative momentum. The R = 1 fm line fails to re- 
produce the data but, as R is increased, the minimum in 
the curve moves towards lower momenta and the data are 
better described. Within the framework of this simplistic 
analysis, the data seem to require a pion- production op- 
erator with a fairly long range component. It was noted 
for the TSL measurements of the differential cross sec- 
tion da I dp [H that the momentum dependence could be 



described by taking the exchanged particle in Eq. (fT2|) 
to have the mass of the pion but that such a description 
was no longer possible when the p-meson mass was used. 



VIII. SUMMARY 

We have measured the angular and momentum depen- 
dence of the cross section and analyzing power of the 
Jlp— >ppir° reaction at an incident energy of 390 MeV. Re- 
coil protons stopped in a scintillation counter have been 
identified by the AE-E technique. Though the counters 
in the horizontal plane covered angles between 15° to 35° 
to the left and right of the beam, the vertical acceptance 
was quite small and this feature is crucial in any the- 
oretical description of the data. Only high relative mo- 
mentum p > 160 MeV/c were registered and, under these 
conditions, final P-state proton pairs play an important 
role. 

Although the dependence of the cross section on the 
pion angle shows a larger anisotropy than in the case 
of data obtained with almost full acceptance @, 0, 0]j 
this is in large part the result of our coplanar geometry. 
After taking this into account by using the IUCF pa- 
rameters @, our results do not contradict those recently 
published by COSY-TOF Q. 

The shape of <ToA y (8 q ) follows quite closely a smd q 
form, which is consistent with this observable being dom- 
inated by PpPs interference. The contribution from the 
SsSd term is negligible at the high relative momentum 
studied in this experiment. However, the behavior of A y 
at small p could be a useful probe to investigate the role 
of the SsSd contribution [28| . 

The momentum dependence of the analyzing power at 
large p is also consistent with the PpPs interference in- 
terpretation. The variation can be explained within a 
very simple model by taking the pion-production opera- 
tor to have a rather long range. 

Our results, taken in conjunction with the IUCF 
double-polarized data Q and other published results on 
the unpolarized cross sections, may provide the extra in- 
formation necessary for the understanding of the pro- 
duction mechanisms for P-state protons. However, to 
succeed in this, a much more sophisticated theoretical 
model is required to replace the rather qualitative ap- 
proach used here to describe the data. 
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